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In dielectric waveguides operating at optical frequencies, the primary 
cause of time dispersion of narrow pulses can be mode conversion. In 
this paper we argue that under certain assumptions a dielectric waveguide 
acts as a linear system in intensity. That is, given the intensity input, the 
intensity output is equal to the input convolved with an intensity impulse 
response. We show that contrary to intuition, the width of the impulse 
response gets narrower when coupling between guided modes increases. 
Using the perturbation results of D. Marcuse, we obtain an interesting 
model of energy propagation down imperfect guides. We conclude that the 
intensity response width increases as the square root of the guide length for 
sufficiently long guides and approaches a gaussian shape for sufficiently 
long guides. 

We conclude from the theory that the dispersion in dielectric wave- 
guides may be orders of magnitude below that which was previously ex- 
pected in guides of sufficiently long length having properly controlled large 
amounts of mode conversion. These theoretical results have not yet been 
verified experimentally. 

I. INTRODUCTION 

In multimode dielectric waveguides operating at optical frequencies, 
the primary cause of time dispersion of narrow pulses can be mode 
conversion. In a geometrically perfect guide with more than a single 
mode, energy initially launched in a given mode remains in that mode 
as it propagates down the guide. Physical guides have imperfections 
from perfect geometric shape (e.g., roughness at the core-cladding in- 
terface of a nominally right circular cylindrical guide) which allows 
energy to couple between modes during propagation down the guide". 
Since group velocities differ in general amongst the modes, a pulse of 
energy initially launched in a single mode or combination of modes 
will be broadened due to the spread of propagation times of different 
parts of the energy. 
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In this paper we argue that under certain assumptions, a dielectric 
waveguide acts as a linear system in intensity as well as in voltage. 
That is, we show that the relationship between the input intensity and 
output intensity of the guide is defined in terms of an intensity impulse 
response. We argue that this intensity impulse response, for sufficiently 
long guides, has a mean-square width about its mean which increases 
only linearly with length. Further, in the limit of very long guides, we 
argue that the response shape is gaussian. We also show that the 
greater the coupling between modes, the less the time dispersion — a 
result which at first contradicts intuition. Finally, we obtain quantita- 
tive results and an interesting model of an optical guide under certain 
assumptions. 

We conclude from the theory that the dispersion in dielectric wave- 
guides may be orders of magnitude below that which was previously 
expected in guides of sufficiently long length having properly controlled 
large amounts of mode conversion. These theoretical results have 
not yet been verified experimentally. 

II. AN OUTLINE OP THE ARGUMENTS 

We next outline the steps of the derivations to follow, so that the 
reader can follow the train of thought. 

We start with the fact that the optical guide is a linear system in 
voltage. That is, if we expand the input signal in spatial modes and 
expand the output signal in the same modes, then the time varying 
coefficients of the modes at the output are related to the coefficients at 
the input by a set of voltage impulse responses. We then make an as- 
sumption about the associated set of transfer functions (Fourier trans- 
forms of the impulse responses) which allows us to argue that the set 
of average output intensities and the set of input intensities are also 
related by a set of impulse responses. Thus the guide is also linear in 
intensity under the assumptions. 

We next argue that for sufficiently long guides, these intensity im- 
pulse responses coupling a chosen input mode coefficient and a chosen 
output mode coefficient are indifferent to the modes chosen except per- 
haps for a magnitude scale factor. 

Finally, this allows us to show that for guides longer than the above 
scale, the intensity impulse response which is now in common for all 
input-output pairs has a mean-square deviation about its mean which 
grows linearly in length, and which approaches the gaussian shape in 
the limit of long guides. 



DIELECTRIC WAVEGUIDE 845 

III. THE OPTICAL GUIDE AS A LINEAR SYSTEM IN INTENSITY 

3.1 The Random Channel 

We now argue that under a simple assumption, the expected value 
of the intensity at the output of a random channel is related to the 
intensity of the input to the channel by a simple convolution with 
an intensity impulse response. We start with input baseband signal 
a(t). We use a(t) to linearly modulate a carrier m(t) which may 
be coherent or a stationary (wide sense) random process centered at 
frequency /„. We assume that the result x(t) = a(t)m(t), has band- 
width B, i.e., its spectrum extends from — B/2 + /o to B/2 + f . 

We pass x(t) through a time invariant filter with a random impulse 
response h(t) representing the channel, resulting in the final output 
y{t). We define a Fourier transform relationship between the function 
A(/) and the function \(t) 

A(/) = f exp [i2TJ(t)]\(t) dt, (1) 

A(/)^X(0- 
By simple linear system theory if 

X(j)^x(t), 

H(j)^h(t), (2) 

Y(j) <=> y(t), 
then 

Y(f) = X(j)H(j). 

Define the envelopes of x(t) and y(t) by 

x(t) = V2 Re ix t (i) exp (*2WoO}« (3) 

y(t) = V2Re {t/.(0exp(i2T/ o «)}. 
The intensity of the input and output signals are defined as 

/ in (0 = wor = « a (0 \™.W, (4) 

/o„ t (0 = |y.(*)l 2 , 

where m e (t) = carrier envelope. 

Assumption: Stationarity of channel transfer function 

(H*(a)H(f + a)) = T(j) (5) 

provided /« - B/2 < «, / + a < /„ + B/2. 
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The assumption, while apparently arbitrary, is essential to the results 
which follow. Perturbation results of Marcuse, 1 to be discussed in 
Section 4.1, indicate that equation (5) may be satisfied for the input- 
output temporal transfer function of a given spatial eigenmode of an 
optical dielectric waveguide with mechanical imperfections, provided 
the mechanical imperfections satisfy constraints also to be discussed. 

Define for any function U(a) 



UM = U(a) a^O, 

= a < 0. 

Then it has been shown that (See Appendix C) 



(6) 



|y.(01 2 «2 f Y+ (f + a)Y*(a) da. (7) 

Then clearly 

Ij/XOI' <=> 2 r X + (f + a)H + (j + a)H*(a)X*(a) da. (8) 

J -co 

Using equation (5) we obtain 

<)y.(0 1 2 ) <=> 2r(/) £ <X + (/ + a)XJ(«)) da, (g) 

<l2/.(0l 2 )«r(/)/ in (/), 

where 

/ in (/) « (T te (0>. 
Thus* 

</.««))= </i»(*)>*7«) (10) 

where 

7 (i) «=* r(/). 

Thus we have a linear system relationship between the channel input 
and output intensities. 

3.2 Extension to Vector Channels 

Suppose we have a vector channel (corresponding to multimode 
guide) consisting of a vector of L input functions 



* The notation x(t) * y(t) signifies convolution: 

/•CO 

x(t) * y(t) — I x(t — u)y{u)du. 

V — CO 



and L output functions 
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The input and output functions are related by an L X L matrix of 
impulse responses 



where we have 



6(0 = [h fi (t)] 

L 

vi(t) = E MO **,(*), 



(ii) 



(12) 



i.e.. 



Y(«) = 6(0 * X(0- 

Now consider a cascade of two vector channels having impulse response 
matrices '8(0 and 2 6(0- The input passes first through channel 1 and 
then through channel 2. The output of channel 2 is given by 

Y(o = 2 e(0 * l e(0 * x(0. (13) 

Define the envelope of y k {t) , Vke(t) : we know that 
\y k .(t)\ 2 <=> 2 f Y t+ (a + j)Yf + (a) da, 

•'-oo 

.(01" ^aflSIZ %!+</ + ")%!+(/ + a)X i + (j + a) 

J -to I / m n 



\Vk 



(14) 



Assumption a. Stationarity of Mode Transfer function, 
b. Mode Transfer functions uncorrected. 

( l H ti {j + a) 2 H kl (j + a yH*M 2 H* n (od) 

= l T ti (J) a r kt (j) 5,.„ 5,, m 6 k . r (15) 
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where 

&x,y = Kronecka delta 

for {/ + a, a) e input signal bands. 

We are implying that randomness, especially in phase, erases corre- 
lation between the transfer functions of different modes. The validity 
of this assumption for optical guides with more than one mode will be 
discussed in Section 4.1. 

It then follows that 

(\yM\ 2 ) <=> 2 f" £ £ l TM 9 T u <fM? t +<j + a)XfM)) da, 

J-« ,-. «-, (16) 

<l»*.(or>« £ E"r„(0 l r„0K w «) 

.-1 1=1 

where 

/wo) ♦♦ <i mo r>. 

Thus under assumption (15) we have 

<]3/*.(oi 2 > = £ £ wo * wo * <Mor> d7) 

j-1 1=1 

where 

7«(«) «=> r„(/). 

Forming the matrix *(?(£) with elements ^kiit) and similarly 2 G(t); 
the vectors of input and output intensities are related by 

| F, | 2 = 2 G(<) * l G(<) * |X. | 2 . (18) 

Thus the vector channel is a vector linear system in intensity as well 
as voltage [compare equation (18) to equation (13)]. 

3.3 A Limit Theorem for a Cascade of Vector Channels 

Now consider a cascade of a large number M of vector channels, 
each behaving as described in Sections 3.1 and 3.2, i.e., if |y r (£)| 2 is 
the vector of the average intensity responses at the output, \X e (t) |- the 
input intensity vector; we have 

\YM* = (£ G >(o) * |X.(0P = Gr * |x.(Q|" 
where 
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U 

G T = * (?,(0 = G u * G^ • • • * G, . (20) 

We would like to argue now that for sufficiently large M, all the 
elements of G T are identical in waveform, differing at most by a con- 
stant. In other words, we would like to argue that the shape of the 
intensity response between any input mode and any output mode is 
indifferent to the choices of input and output modes, for sufficiently 
long guides, except perhaps for the magnitude of the responses. 

We shall prove our results for a lossless two-mode guide. Define 
yij(t, b) as the j -> i intensity impulse response for b sections of 
guide. We obtain (see Appendix A for derivation) 

?«(*, &) - S Tuft 6 - «) * Tm(0 * T 2 * 2 B_1 (0 (21) 

R = l 

where 

y&*-\t) 4 y 22 (t) * y 22 (t) ■ ■ ■ R - 1 times 
and 

y„(t 0) = 5(0- 
Now assume that the guide is lossless, i.e., 



/: 



m lTu(0+7ii(0l* - 1- (22) 

Further define 

7 .-,-(«) - aapiiit), 
r p „{t)dt=l, 0<o„<l. (23) 

J -00 

Thus 

7ii(*. b) - S Yu(«, 6 - «) * p»i(0 * P?»*" I (0(a.ia.V 1 ). (24) 

R-l 

For the lossless guide, and b sufficiently large -y n (6 — R) « 7n& &) 
for # « b. Furthermore a convolution of 7„(i, b — R) with 
p»i(0 * P*2 R ~'(0 is approximately equal to 7n($, b — R) for /? « 6 
since the response 7n(£, & ~" ^)> which is a convolution of b — R terms, 
has a narrow spectrum compared to the other R term convolution for 
6 » R. Furthermore a^a*^ 1 — > for R large. Thus for b sufficiently large 

7».(0) «7iiG,fr)(««/(l - a 22 )). (25) 
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Similarly we have 

Tuft b) - E-Yiaft b-R)* pM * rfi^'CO*!*", 



&y 12 (t,b) T ^— = y l2 (t,b), (26) 

1 On 



(since a 2 i + a„ = 1). Thus 

_1 1 

Oft b) S p («, 6) 



1 + 7, 1 + 7, 

1 1 



= P {t, b)k, (27) 



.1 + 1/7, 1 + l/,_ 

where n = a 21 /a l2 , and p(£, 6) = y n (t, &)//"„ 7nft &) <8- Finally we 
obtain the response of a guide of kb sections 

Gft fc&) = P **(*, 6)A. 

Note that A is idempotent, i.e., A 2 = A and P {t, b) is a positive unit 
area function. 

3.4 Application to Long Optical Guide 

For a multimode lossless of guide sufficiently long length, I = kL, 
with finite coupling between all modes, we can generalize equation 
(27) to conclude that the intensity impulse response between an input 
and output mode is a constant times some positive unit area function 
p(t, L) convolved with itself l/L times when L is a scale on which 
equation (27) holds in the generalized case (more than two modes). 
Since the central limit theorem states that the convolution of a large 
number of unit area positive functions approaches a gaussian shape,* 
we conclude that the impulse response should approach a gaussian 
shape in the limit of long guides. Further, the impulse response's 
second moment about its mean increases linearly with increasing guide 
length for guides longer than L. f That is, the second moment about 
the mean of the response is 

Af t © = M 2 (L)l/L. (28) 

If T! is the "differential delay" (time/meter) of propagation in the 



* Provided that S" a p(t, L)t 2 dt < oo, when we add similar independent random 
variables with finite second moments, the probability density of the sum, which 
is the convolution of the individual densities, approaches a gaussian shape. 

t The second central moment of a convolution is the sum of the individual 
second central moments. 
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slowest modes and t 2 is the differential delay in the fastest mode, then" 

MiiL) s k^fH m !&£. (29) 



Therefore, 



M 2 {1) < ^- LI 



for any L where equation (27) holds [for an N mode case we have an 
N X N matrix multiplying p(t,l)]. 

IV. QUANTITATIVE RESULTS 

4.1 Perturbation Theory 

We shall now apply the above results to the case of a lossless slab 
dielectric waveguide previously studied by Marcuse. 1 We expand the 
input field to the guide as 

<t, X, 0) = L e k (t, 0)fc(x) (30) 

where x is the cross-sectional position parameter and the \J/ k (x) are 
the eigenmodes of the guide. The field a distance I down the guide is 
written as 

e(t, X, I) - E *ft WkC*)- ( 31 ) 

We have a linear voltage impulse response relationship between the 
vector of input voltages [e k (t, 0)] and the vector of output voltages 
[e k (t, 1)]. Defining E k (ta, I) as the Fourier transform of e k (t, I) we have 

E h (v, = £C„(<-,9#,(«,0). (32) 

Marcuse has shown that a perturbation theory solution for the C ki (w, I) 
is given by 

C ki (o, I) = Ki exp [ifi t (p)t\ [ g{z) exp [i[fi t (u) - /3y(«)» rfz (33) 

Jo 

where X fc , is a constant weakly dependent upon co and g(z) is the wall 
perturbation from straightness. It is assumed k ^ j [For fc = j, 
C tj («, = 1]. We have therefore 

* The right side of equation (29) is the mean-square intensity impulse response 
width if the response consists of an impulse of area i at the shortest delay and an 
impulse of area I at the longest delay. 
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C„(«, OC* (CO + <7, I) 

= X 4 ,X* exp [i(Pi(u) - |8,(« + a))l) f J g(z)g{z') 

•exp [i[(p k (w) - £,(«))* ~ (&.(« + «) ~ fip(f* + a))z']) dz dz' . (34) 
Defining the correlation function 

(Q(z)g{z')) = R,(z - z') (35) 

(we assume g(z) is a wide sense stationary process) we obtain 

<tf w («, Z)C* (« + cr, 0) = *«*& exp (» Aft?) J J R„(z - z') 

•exp [i(fi k (u) - Pj<a)(z - z')] exp [i A/3(co, cr)«'] dz dz' (36) 
where 

A/3 = (ft(«) - B (co + cr)) - 08,(0,) - ft(« + a)), 

A/3 2 = ft(«) - 0> + cr). 

If jR b (z — 2') drops off quickly for (z — z') in an interval of length I, 
then we have the approximate result 

(O, QC7* (« + *, 0) ac \MlSMa) - ft(«)) 

•expfr-^O 1 "^^^ (37) 
i Ap(co, cr;t 

where £„(■) = Fourier Transform of #„(•), and is assumed to be 
constant as a function of /3*(a;) — /3,-(w) for « within the excitation 
bandwidth. For k = n, j = p 

« - &>• 

That is we assume o- is small enough so that there is negligible disper- 
sion of energy travelling in a single mode. Thus, the intensity impulse 
response between input ; and output k is (see Fig. 1) 

y ki (t, I) ~ K,KA(0M - PMW - r f l) (39) 

where 

f(t) - 1, H0,Ar mi l]; 

= 0, otherwise; 
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for I small enough for the perturbation theory to hold. 

From equations (3G) through (38), it is clear that for cases where we do 
not have k = n, j = p the correlation function C t ,(co)C„*(w + a) will be 
negligible provided A/3(a>, a) is sufficiently large in the band of input 
frequencies. Thus we can use the perturbation length impulse response 
to find a long-guide response by means of equation (20). 

From equation (37) we see that we increase coupling between modes 
by making the mechanical perturbation spectrum large at frequencies 
which correspond to the difference of the inverses of the phase velocities 
of the modes at excitation frequencies. It should be emphasized that 
while making the perturbation spectrum high at frequencies that 
couple guided modes, we will wish to avoid making it too high at 
frequencies that couple guided to unguided radiation modes since such 
coupling results in loss. 

4.2 A Hydraulic Model of Dispersion 

We shall now show that the perturbation results imply a model 
which is an interesting interpretation of the propagation process, and 
which allows easy computation of the response of a long guide. 

Suppose energy traveled down the guide as follows. We start with a 
large number of indivisible bundles of energy at the guide input. Each 
bundle begins propagating down the guide randomly jumping from 
mode to mode. At any point down the guide, a bundle travels at the 
group velocity associated with the mode it is currently in. At any posi- 
tion, the probability that a bundle will jump to mode k, given that it 
is in mode /, in the next increment of distance dl is X 4 ,X$S B (0 fc (ci>) — 

Since we have a very large number of bundles, the output response 
of the guide in intensity should have the same shape as the probability 
distribution of the arrival time at the output of an individual bundle. 
For a short guide of length L, the probability that a bundle is in mode k 
given that it started in mode j is | \ ki | a S (p k (u) - /3,(o>))L and its 
arrival time distribution is given exactly by Fig. 1 . Since this distribution 
is the perturbation solution for the intensity impulse response of a 
short guide, we see that the hydraulic model gives the same result as 
the perturbation theory. Further, a little thought will show (see Ap- 
pendix B) that the extrapolation from a short guide to a long guide in 
the hydraulic model is analytically the same as equation (20). Thus 
any technique which can be used to determine the intensity impulse 
response characteristics using the hydraulic model will be valid for the 
solution of equation (20) using the perturbation results. 
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Fig. 1 — Output intensity response for modes ;', m and k, given mode ;' is excited. 



4.3 The Solution of Some Hydraulic Model Response 

4.3.1 Characteristics 

We now wish to determine some probability density moments of 
the arrival time of a bundle of energy at the output of a long guide 
recalling that this has the shape of the guide impulse response. 

Let H{V) be the mode a bundle is in at distance V down the guide. 
In that mode the bundle travels with differential delay (time/meter) 
th = t{V). The total propagation time down the guide is 



(40) 



(41) 



T = f t(1') dl f . 

Jo 

The expected propagation time down the guide is 

(T)„ = f (t(Z') dl')„ . 
The variance about the mean is 
<(T - (T) av )\ v = (£ £ (r(Z') - (r) av (l'))(r(l") - (r)Un) dl' dl"} 

= \J J KM, I") dl' dl"\ - (T)l (42) 

where R T (l', I") = E{t{1')t{1")). We need the correlation function 
R T (l, I') and the mean (7-) BV (0- 
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4.4 Calculation jor a Lossless Two-Mode Guide 

For the lossless two-mode guide we have an energy bundle making a 
Poisson number of mode changes in any length L with mean 
I Xi2 T *S (/3,(w) - /3 2 (w))L. The correlation function R T (l, V) (assuming 
we start off randomly in one of the modes) is that of a random tele- 
graph wave 2 and is given by 

R T (l, I') = \ |Ar 12 | 2 exp (-2 |J - V\/l e ) + i^^j H " 0" (43) 

where 

l/lc = I X.2 | 2 ^(^(o>) - 2 (io)), 
and 

Ati2 = tar " fo) m Ti " T2 • 

We obtain the mean and second moment about the mean of the 
intensity response of a guide of length L. 



<r>„ = {^)l, 



<(T - (T)„)%, = ^ T 'f l - L \l - i (1 - exp (-2LA))] , (44) 

um «r - <r>„) v - ^p!* , 

where 

L = [lx,»l a S.03.(«) -AW)]" 1 . 

Um <(T - (TUV = ^£¥ , 

(compare equation (44) to equation (29)). 

4.5 Extension to the Two-Mode Guide With Loss 

We can use the hydraulic model to extend the above results to a 
two-mode guide with loss and differential loss. 

Assume that when travelling at distance dl in mode ;', a bundle of 
light has probability ajdl of being absorbed. 

If in travelling down a guide of length L, the bundle spends a dis- 
tance Lj in mode 1 and L 2 in mode 2, then the probability that it is 
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not absorbed is 

P = exp [— (ctjji + a 2 L 2 )\ = exp [-((a, - a a )L, + a 2 L)]. (45) 

The number of bundles entering the guide at time zero and arriving 
at the output end at time t in the presence of loss equals the number 
that would arrive at time t in the absence of loss times the probabil- 
ity that a bundle with total travel time t is not absorbed. But we have 

L\ ~r L 2 = L, 

t = T X L X + T 2 L 2 = (r, — T 2 )Lx + T 2 L, 

P = exp [-foL, + a 2 L 2 )] = exp |-(a« -^7^ + « 2 ^)} 



(46) 



where 

At 12 = Tj — t 2 Aq; = a, — a 2 - 
With a little algebra we obtain 

P = exp [-{^T « " (t>«vL) + <a), lv Z,}] (47) 

where 

(«) nv = (a, + a 2 )/2, 
<r> nv - (r, + r 2 )/2. 

Thus the intensity impulse response for a two-mode guide with loss 
is equal to the lossless response multiplied by P of equation (47). 

Note that the gaussian shape for long guides still holds because the 
product of a gaussian and an exponential envelope is a shifted 
gaussian. 

V. CONCLUSIONS 

We can conclude at least one important result. Long optical fiber 
waveguides need not have large dispersion due to random imperfec- 
tions if properly controlled mode coupling exists. From equation (44) 
we see that a mechanical perturbation spectrum which is peaked at 
frequencies that couple guided modes will lower dispersion. However, 
to avoid loss, we must not make the mechanical perturbation spectrum 
too high at frequencies that couple guided and radiating modes. 
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The above conclusions have been obtained by D. T. Young and 
H. E. Rowe, 3 for the two-mode guide by solving the coupled line 
equations directly under the assumption of white noise coupling. 

APPENDIX A 

We wish to derive equation (21) from the following relationship 



G(t, b) = [ 7ii (t, b)} = [ yii (t)]*[y<MY 



(6 Times), 

h 3 = 1, 2. 



(48) 



Equation (48) implies the following model shown in Fig. 2. The 
transfer function y 2 i(t, b) is the overall transmission response between 
input 1 and output 2. This can be obtained by adding up the trans- 
mission responses over all different paths between input 1 and output 2 
using any desired bookkeeping scheme. Every path between input 1 
and output 2 must pass through the y 2i (t) function jor the last time in 
some section. If a path passes through the y»i(t) function for the last 
time in the fifth section from the end, then it must pass through four 
7 22 (0 functions on its way to output 2. The sum of the path transfer 
functions between input 1 and the input to the 721 (0 function in the 
fifth section from the end is y u (t, 6 — 5). Thus the contribution to the 
overall transfer function between input 1 and output 2 due to all paths 
which pass through a 721 (0 function for the last time in the fifth section 
from the end is ?„(*, 6 - 5) * y 2i (t) * y* 22 (t). Equation (21) merely 
expresses the sum of the contributions over all positions of last passage 
through a 721 (0 function. 




•b SIMILAR SECTIONS 



Fig. 2 — b-section guide. 
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Y«(*, b) = £ Ynfo 6 - R) * y 2l (t) * ytf- l (t). (21) 

R-l 
APPENDIX B 

We wish to show that the equations for obtaining the intensity 
response of a long guide, given the intensity response of a short guide, 
for the hydraulic model are identical to the extrapolation equations 
for the intensity response given by equation (20) . 

Suppose we have the probability density, for a short guide, that 
a bundle of energy starting off in mode ; of the guide (at time 
zero) arrives at the output of the guide at time t in mode i. Thus we 
have the matrix of densities P{t, L) where L is the guide length 
and the elements py(t, L) are the previously described densities. 
Let Ij(t, 0) be the probability density that a bundle of light arrives 
at input ; at time t. Let Ij(t, L) be the probability density that a 
bundle arrives at the output position L in mode ; at time t. Let I(t, • ) 
be the corresponding vectors. Using the laws of addition of random 
variables we obtain 

W,L) = T,pUt,L)*I.(t,0) 
of in matrix notation 

I(t, L) = P(t, L) * I(t, 0) 
therefore 

1(1, kL) = (*P(t,L))*I(t,0). 

We see that the probability density of the output arrival mode and 
time of a bundle of energy for a long guide, which corresponds to the 
intensity response, is extrapolated from the short-guide response 
exactly as in equation (20). Thus since the perturbation results of 
Marcuse correspond to the hydraulic model in the limit of short guides 
and satisfy the conditions for extrapolation using equation (20), it 
follows that properties of the hydaulic model solution for long guides 
will correspond to the solution of equation (20) starting with these 
perturbation results. This is true no matter what techniques we use 
to find these hydaulic model properties. 

APPENDIX C 

We wish to establish that for a narrowband high frequency signal 
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y(t) = (y.(t) exp (W) + y*(t) exp (-*W))/V2 

of carrier frequency / = co /2ir and envelope y.(t), the intensity is 
given by 

\y.(fi\ % = 2 / [j"° FJ(f)F*tf + a) dfj exp (-t2«ri) da, 

t=>2 f Y*(j)Y + (f + a)dj. 
Define 

r.(fl = f VM exp (t2ir/<) dt, = v2 7 + (/ + /„), 

J -oo 

[provided y,(£) is narrowband compared with / ] ; 

\yXt)\ 2 = f° Y.(D exp (-2iw1t)Y*(r) exp (Oaf Q d/ df, 

J -to 

= f f exp 1-iMf ~ r)t][Y.(1)Y*(n) dfdf, 

= f° f exp {-i2mt)[Y*{r)YXV + 7)] dy df, 
where y = / — /'; 

- 2 f" exp (-*2ryi)[Fl:^)y + (ff + 7) dy dg] 

J -co 

where g = (/' +/ ) . 

Q.E.D. 
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